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k*" ! Abstract 

t — ■ We will study open and hidden charm scalar meson resonances within two different 

models. The first one is a direct application of a chiral Lagrangian already used to 
study flavor symmetry breaking in Skyrme models. In another approach to the 
^sO • problem a SU (4) symmetric Lagrangian is built and the symmetry is broken down 

to SU{3) by identifying currents where heavy mesons are exchanged and suppressing 
those. Unitarization in couple channels leads to dynamical generation of resonances 
in both models, in particular a new hidden charm resonance with mass 3.7 GeV is 
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Oh. predicted. The small differences between these models and with previous works will 

be discussed. 



X 1 Introduction 
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The discovery and soon after confirmation of charmed scalar resonances by BaBar 
and Belle pQ, [2], [3] has opened a controversy about their structures. In the qq 
picture these resonances are naturally assigned as 3 Po states in the 2S+l lj notation, 
but calculations done long before in the framework of quark model potentials [4] had 
mass predictions which turned out to be more than 100 MeV off the real mass of the 
states. Lattice calculations also fail in calculating the masses with a qq assignment 

m- 

This situation has sparked the discussion whether these resonances could have 
a different structure. Some authors have suggested a qqqq structure [6], [7] or a 
mixing between four quarks and the usual qq structure [8]. Also molecular states 
have been suggested [9], [10] . [TT] . For discussions on these and other controversial 
heavy mesons refer to [12] . [13j . [14] . 

Unitarized coupled channel models have also been considered for the study of 
these resonances in [15] . |16| and [17j . These works have used a chiral Lagrangian 
based on heavy quark symmetry [18] . [19] . [20] . |21j for the open charm sector which 
neglects exchanges of heavy vector mesons in the implicit Weinberg- Tomozawa term. 
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We intend to extend the study for all possible sectors of the interaction, including the 
hidden charm and the double charmed sector. The exchange of heavy vector mesons 
is also taken into account in our approach, with the corresponding terms properly 
accounting for the larger mass of the heavy vector mesons. 

In the present work we will construct a Lagrangian for the interaction of the 15- 
plet of pseudoscalar mesons in S I J7(4). SU(4) symmetry breaking will be considered 
by suppressing the exchange of heavy vector mesons. The SU{3) structure of the 
interaction will be thoroughly analysed and unitarization in coupled channels will lead 
to the generation of scalar resonances corresponding to poles in the T-matrix. For 
comparison, in the open charm sector, we will also solve the problem with yet another 
Lagrangian which has been considered in the study of flavor symmetry breaking 
effects in Skyrme models. This Lagrangian gives similar results, supporting our 
present model. 

The work is organised as follows: in the next section a brief review of the structure 
of the SU (4) 15-plet will be presented. Section 3 is dedicated to the explanation of 
the construction of the Lagrangian and in section 4 the theoretical framework for 
solving the scattering equations in a unitarized approach is presented. Section 5 is 
dedicated to analysing the results and a brief summary is presented in section 6. 

2 The 15-plet 

In this work the framework already used to study the interaction of the octet of 
pseudoscalar mesons in SU(3) [22] will be extended to include charmed mesons. 
This will involve some extrapolation to SU(4). In the qq picture, mesons involving 
charm will be classified as 4 ® 4 = 15 © 1, hence belonging to a 15-plet or a singlet. 
It is interesting to understand how the 15 representation of SU(4) breaks down into 
representations of SU(3) and in which channels the interaction of the multiplets will 
be attractive or repulsive. 

The pseudoscalar mesons are represented by a 15-plet of SU(4) as shown in figure 
HJ Once SU(4) symmetry is broken into SU(3), the 15-plet breaks down into four 
multiplets of the lower symmetry, an octet, a triplet, an antitriplet and a singlet: 

15 — >1©3©3©8. (1) 

The octet and the singlet have null charm quantum number, the triplet and the 
antitriplet have negative and positive charm quantum number, respectively. When 
studying the meson- meson interaction, one can decompose the scattering of two 15- 
plets of SU (4) according to its SU (3) inner structure, see Table [Q 

The interaction 8 © 8 is already very well studied [22] -[29]. It is already known 
that in s-wave it generates dynamically poles in the T-matrix which are associated 
with the /o(980), k, ao(980) and a resonances. Also the interaction 3® 8 has already 
been studied [15] . [16] . |17j in some different approaches. In this sector, as we will 
show in the following sections, the interaction, when diagonalized in a SU(3) basis, 
is attractive in 3 and 6, while repulsive in the 15, so one could expect to generate 5 
poles in (S, 1): from the 3 with isospin and \ and strangeness 1 and 0, respectively; 
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Figure 1: 15-plet from 57/(4) with its mesons assignments. 



charm 


Interacting multiplets 


2 


3(8)3 -> 3©6 


1 


3©8^15©3©6 
3© 1 -> 3 





3©3 -f 8© 1 
1® 1 -> 1 
8© 1 -> 8 

8 ® 8 -»• 1 © 8 S © 8 a © 10 © 10 © 27 



Table 1: SU(3) decomposition of the meson- meson interaction in 57/(4). The sectors not 
shown in the table correspond to the C = — 1, —2 states which are just charge conjugate 
states (antiparticles) from the ones shown. 
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from the 6 with I = 0,^,1 and S = —1,0, 1, respectively. Moreover the interaction 
in the 3(8)3 is attractive in both the 8 and the 1, so one can, in principle, expect four 
new resonances in the hidden-charm sector. The interaction in the C = 2 sector is 
repulsive in the 6 and in the 3 the interaction vanishes. 

Apart from studying the different sectors separately it is interesting to see how 
the mixing of states from different sectors with the same SU(3) representation of 
Table [U affects the interaction. 

Furthermore, if heavy resonances are generated from the 3<8> 3 one can expect, in 
principle, that the mixing of those heavy channels with light ones coming from 8 (8> 8 
will make its width quite large because of the large phase-space available for decay. 
However, we shall also see that there are subtleties in the interaction which suppress 
these decays. 



3 The Lagrangians 

The SU(3) lowest order chiral Lagrangian reads [30], [31] : 

f 2 f2 2 

C x = ^Tr (dpUdPXJ) + J -^-Tr (u + tf- 2) (2) 

where U is the field containing the pseudoscalar mesons from the SU(3) octet and 
Tr represents a trace in flavor space.: 



U 



e fa 

V2^ V6 



V 



7T 

K- 



7T+ K+ \ 

-tt° 1 n v-o 



(3) 
(4) 



Flavor symmetry breaking effects can be introduced with two new terms in the 
Lagrangian [32], [33] : 



C 



SB 



I, 



f 2 K m K ~ fl m l 



6 

fx ~ fn 

12 



Tr((i-V3\ 8 )(U + U^-2)) + 



Tr((i- v / 3A 8 )(C/y^ + y' 1 [/ t ) 



(5) 
(6) 



where As is the SU (3) generator. 

In [34] these Lagrangians are extended to SU(n). In this new approach the 
symmetry breaking sector is written as: 



fc=3 



~ J2 hTr 1(1 -J \k(k - l)A fc2 _ 1 )(C7 + - 2) 



(8) 
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but now U belongs to a SU(n) representation. 

By expanding the U matrix until fourth order in the meson fields, one can identify 
the mass and kinetic terms for each field and fix the symmetry breaking parameters 
for S*7(4) and SU(3) as: 



73 = 




(9) 


h = 




(10) 


74 = 


\{fl + f 2 K-Vl) 


(11) 


<5 4 = 


j-2 2 1 f2 2 1,2 2 
/#D - ^K m K ~ ^L m n 


(12) 



In this work we will consider only the difference between fjj and f n which is about 
70% and we will make the aproximation fK = fn- 

For constructing our model we will first consider a SU(4) field containing all fields 
from the 15-ple!0: 



( _i_ JL _i_ Jk_ 

V2 VE Vvi 



7T 

D° 



TT 1 

-7T° I JL _|_ JISL. 

K° 

D+ 



K+ 
K° 



Df 



D° \ 
D~ 



(13) 
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A current is then defined: 



(14) 



and a Lagrangian is build by connecting two currents and adding an extra term 
proportional to the square mass of the fields: 



Cpppp = ^Tr{J^ + $ 4 M). (15) 

SU(4) and SU(3) flavor symmetry breaking already arise from the mass term when 
the matrix M is not proportional to the identity matrix. We take: 

^^What here is called r\ and r\ c are actually 77s and 7715 states, which mix with a singlet, r\\ to form 
the physical states rj, rj' and T) c , but in this work this mixing won't be taken into account, and it will be 
considered that the physical states are just described by their most important components. 
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2m 2 D - 


ml ) 



(16) 



The term with the matrix M is exactly the same one appearing in the chiral 
Lagrangian of eq. after breaking SU(4) and SU(3) by means of ([5]). The term 
To]) appears for four meson fields from the kinetic term, dJJd^U , of 



J M J M in eq. 

the chiral Lagrangian in eq. ([2]) if U is taken as a SU(4) representation by means of 
replacing (/>$ in eq. (j4]) by in eq. (fl~3l) 

We will also implement a different sort of symmetry breaking in a way that we 
explain below. 

The constant / appearing in the Lagrangian (|15p is, in principle, the pion decay 
constant (in this work f n = 93.0 MeV). However, a different one will be used for 
the heavy mesons. In this latter case, the f 2 appearing in the amplitudes should be 
thought as the product of y/J for each meson leg in the corresponding vertex, with 
/ = f n = 93.0 MeV for light mesons and f = /d = 165 MeV for heavy ones. This 
value for fjj is of the order of magnitude expected from the experimental point of 
view |36| and lattice calculations |37| . 

Directly applying Feynman rules to obtain transition amplitudes from this La- 
grangian would be too much of a simplification. Indeed, the term J^J^ of the chiral 
Lagrangian is usually visualized as the exchange of a vector meson between pairs of 
pseudoscalar fields in the limit of q 2 « my (the Weinberg Tomozawa term). In 



this case the kinetic term of the Lagrangian of eq. (|15p is SZ7(4) flavor symmetric 
and therefore implicitly assumes equal my for all the exchanges of heavy and light 
vector mesons. In refs [15) . [16] . [IT] an SU(3) version of the interaction based on 
heavy quark symmetry is used which would correspond to allowing the exchange of 
only light vector mesons in the Weinberg- Tomozawa Lagrangian described by the 
derivative term of eq. (|15p . and neglecting the M term [35]. In the present work we 
shall go one step further by allowing also the exchange of heavy vector mesons but 
weighted by their respective squared masses and we shall also keep the mass term as 
done in [32] , [33] , [34] • In order to implement this we first decompose the $ field into 
its SU(3) components: 



(17) 



The I3 is the 3x3 identity matrix and the fields <pi contain the meson fields for 
each i-plet of SU(3) into which the 15-plet of SU(4) decomposes: 



/ 4 + JL 
K~ 



7T+ K+ \ 
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fa = ( D° D+ Df ) 

01 = Vc 

In this way the Lagrangian in (|15p can be decomposed into six parts: 



£>PPPP = ^2/2 + -^3 + -^31 + £83 + £831 + £mass) (18) 

C 8 = Tr(J ss ,4 8 ) (19) 

^3 = J- 33 ^ 3 + Tr(j 3 - 3 ^ (20) 

^31 = 7^31^13 ( 21 ) 

£ 83 = 2(jg 8 ^ + rr(jgg^)) (22) 

^831 = — ^ (^3lAt J 83 + J 38m J 13) ( 23 ) 

£ m ass = 7V(m$ 4 ) (24) 
where the currents are defined as J^- = (d^(j)i)4>j — tfi^d^cfrj). 



Now the exchange of charmed (heavy) vector mesons can be easily identified in the 
different pieces of the Lagrangian by identifying currents carrying explicitly charm 
quantum number. The term accounts for the exchange of light vector mesons 
only. In £§3 the first term is mediated by heavy vector mesons and the second term 
by light ones, £s3i an d £31 have only contributions from heavy vector mesons and 
£3 will still have to be worked out further. 

The separation of the heavy vector meson contribution from £3 is more subtle 
because the exchange of a heavy hidden charm meson in this sector occurs in charge 
and flavor conserving hadronic currents, where the po an d o> also contribute. The 
strategy followed here is to construct a Lagrangian connecting the current with a 
vector field W [35], [38]: 



C PPV = --^=Tr([d^]Vn. (25) 

Here is a 4x4 matrix with the same structure as but with the 15-plet of vector 
mesons instead. The heavy vector meson which can be exchanged in charge and 
flavor conserving hadronic currents is the J^. 

The J^p contribution can be calculated from the Lagrangian ([25]) and it is easy to 
see that when the vector mesons are connecting equal hadronic currents one has a 
contribution with weights 3 and | for light vector mesons and the J^, respectively, 
while the weights are — | and | in terms connecting different currents. Appendix C 
shows in detail how to work out £3. 
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With all these considerations the full Lagrangian can now be rewritten in terms 
of the correction parameters: 




(26) 

'-■> - + ^1 — I (27) 



^ = ~ + — \ • (28) 



Here mi, and run are masses of light and heavy vector mesons. They will be set 
to rriL = 800 MeV and ran = 2000 MeV. With these ingredients the full corrected 
Lagrangian can be written as: 




Note that from eq. (j29[) we can recover the usual lowest order chiral Lagrangian 
for SU(3), which is the term proportional to Tr{ J%%^ Jgg), while the Lagrangian used 
by Kolomeitsev [15] and Guo [17], based on heavy quark symmetry [18], [19] . [20] . 
[21j . is proportional to the term Tr(J^ 3fl J 88 ). Our model has also terms for the 
interaction of heavy mesons only, proportional to rr(J 3 g J^) an d J33 Jg 3 and all 
the other terms are corrections that can be controlled by the parameter 7. 

From this Lagrangian, applying the usual Feynman rules, the transition ampli- 
tudes in appendix A are calculated and used as potential for each possible reaction. 
This potential will be used as the kernel for solving the scattering equation. 

In order to support our results also the chiral Lagrangian with the flavor symmetry 
breaking pieces will be used to solve the scattering problem in the open charm sector. 
Very similar results are found and will be discussed in section 5. 



4 The Scattering Problem 

The amplitudes needed, M(s,9), are written in Appendix A for the Lagrangian of 
eq. (|29p . Since we are only interested in s-wave meson- meson scattering, we first 
project the amplitudes over s-waves, by making a simple angular integration. After 
projecting the amplitudes for s-wave they will be transformed to isospin basis and 
inserted into the Bethe-Salpeter equation which in the on-shell formalism of [22], [39] 
is reduced to an algebraic equation: 

T = V + VGT. (30) 

In this equation V is the potential, a matrix constructed with the tree level transition 
amplitudes for each one of the possible channels, projected over s-wave. The matrix 
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G is diagonal with each one of its non-zero elements given by the loop function for 
the two particles in each channel: 



,2 



V r ill, ;/t 9 — i/L-,-r^ T ui, 9 

on + Log^ + 2 , 1 Log^ + 

JL ( Log !Z^g4 + Log J±H^+^L) ) . (32) 

P in equation (|31|) is the total four-momentum of the two mesons in channel i and 
mi and mi are the masses of the two mesons in this channel. The expression in 
eq. (f32|) is calculated using dimensional regularisation. Over the real axis p is the 
three-momentum of the mesons in the center of mass frame: 



yj(s - (mi + m 2 ) 2 )(s - (mi - m 2 ) 2 ) , . 

P = 2jS • (33) 

In the complex plane the momentum p is calculated using the same expression. 
Eq. (|30p with eqs. (|31H32p makes implicit use of dispersion relations in which only 
the right hand (physical) cut is considered. It was proved in [30] that the left hand 
cut provides a moderate contribution, and more important, very weakly energy de- 
pendent, such that its contribution can be easily accommodated in terms of the 
subtraction constant that we use, in the range of energies of interest to us. 

In this work we will set the loop parameter in eq. f|32|) to fi = 1500 MeV and fit 
the substraction constant, free parameter. 

This loop function has the right imaginary part to ensure the unitarity of the 
T-matrix |23j : 

I m ( Gil ) = -^=. (34) 
Equation (|30p can be easily inverted: 

T = (1 - VG^V. (35) 

When looking for poles in the complex plane one should be careful because of 
the cuts of the loop function beyond each threshold. Bound states appear as poles 
over the real axis and below threshold in the first Riemann sheet. Resonances show 
themselves as poles above threshold and in the second Riemann sheet of the channels 
which are open. 

Over the real axis the discontinuity of the loop function is known to be two 
times its imaginary part [41] so, knowing the value of the imaginary part of the loop 
function over the axis, eq. (|34p . one can do a proper analytic continuation of it for 
the whole complex plane: 

Gil = G i + i^, lro(p)>0. (36) 
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lm(G) 





Figure 2: Upper left: Imaginary part of the loop on the first and second Riemann sheets 
superposed. Upper Right: Real part of the loop in the first Riemann sheet. Botton left, 
right are the imaginary part of the loop in the first and second Riemann sheets, respectively 



G 11 and G 1 refer to the loop function in the second and first Riemanian sheets, 
respectively. 

Figure [2] shows some plots of the loop function in the complex plane. 



5 Results 

The amplitudes listed in appendix A are in a charge basis. First, they are transformed 
to an isospin basis and then to a SU (3) basis by means of the isospin and SU (3) states 
given in appendix B. The SU{3) symmetry breaking is then studied. The amplitudes 
in a SU(3) basis show in which sectors the interaction is attractive (where it may 
generate resonances). The most important term in each amplitude is the s term, 
so when an amplitude has a negative factor multiplying s it is considered to be 
attractive. The results of the diagonalization of the interaction in SU (3) basis are as 
follows: 

8 8 -> 1 © 8 a 8 10 10 27: The interaction here is repulsive in the 
27; there is no interaction in the 10 nor in the 10; in one of the octets, because of 
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its symmetry properties (here we should consider scattering of identical particles if 
SU(3) symmetry is restored) there is no interaction in even I partial waves; in the 
other octet and in the singlet there is attraction, which will lead to the formation of 
4 states, to be identified as the light scalar resonances, a, fo, ao and k. 

3(g) 3 — > 3 ©6: Here there is no interaction in the 3 when the correction factors are 
set to 1 , otherwise the interaction has a p-wave structure (t — u). In the sextet the 
interaction is repulsive, therefore, no double charmed scalar resonances are expected 
from our model. 

3 © 3 — > 8 © 1: The interaction is attractive in both the octet and the singlet 
if the correction parameters (^3, 1^5) are set to 1. In this case, where the large 
mass of the is disregarded, one can see resonances generated. However, since the 
terms with the heavy vector meson have the largest weight in the amplitude, when 
the correction parameters are considered to take into account the different masses of 
the exchanged vector mesons, the resonances disappear for the octet. The singlet is 
always attractive irrespective to the correction parameters. 

3 © 8 — > 3 © 6 © 15: In the anti-triplet and sextet there is attraction while in the 
15-plet there is repulsion. We generate in our model five resonances with charmed 
quantum number, two from the antitriplet and three from the sextet. 

We discuss below the free parameters of the theoretical framework and how we fit 
them. The parameters changed are an and cxl- The a's are the subtraction constants 
for the loop functions. The parameter an was chosen for channels involving at least 
one heavy pseudoscalar meson and a different one for channels where there are just 
light ones, «£. 

One of the novel aspects of the present work is that we allow the mixing of the 
light mesons with the heavy ones in the search of zero charm or hidden charm scalar 
mesons. The first interesting result is that the influence of the heavy meson sector in 
the generation of the light scalar resonances (a, fo, ao, k) is negligible. For instance 
it was checked that different values an have very small effect over the pole position 
for the light resonances. Varying an between -0.3 and -2.3 has less than 10% effect 
over the pole position of the fo resonance, for example. So the heavy sector can be 
worked independently from the light one. 

With this in mind the open-charm (C=l) sector was used to fit an so that the 
position of the pole in the S=l, 1=0 sector match the D* (2317), which has already 
been suggested as being dynamically generated in [15] , [16] and [T7] . After fixing the 
heavy parameter, the ai was fitted by locating the pole position in the sector C=0, 
S=0, 1=1, which correspond to the ao resonance. We also made the fit of an for the 
model involving the chiral Lagrangian. The results are as follows: 

Phenomenological model: a#=-1.3 and ol=-1.3 

Chiral model: a#=-1.15 (we only applied this model for the open charm sector.) 

5.1 SU(3) Symmetry Breaking 

In our phenomenological model it is assumed that the SU(3) flavor symmetry break- 
ing arises from the different masses of the interacting mesons. The mass used for 
each member of the 15-plet is: 
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777^=138.0 MeV, 777^=495.0 MeV, m,,=548.0 MeV, 

m D =1865.0 MeV, 777^=1968.0 MeV and 777^=2979.0 MeV 

Note that there is no isospin breaking in the model, all particles in a same isospin 
multiplet are considered to have the same mass. So, in this work, the Bethe-Salpeter 
equation, eq. (f30|) . was solved with V in isospin basis. 

SU(3) symmetry can then be gradually broken by means of a symmetry breaking 
parameter x which takes values between 1 and 0, 1 meaning symmetry broken as we 
see it in the real world, and symmetry restored. The masses of the mesons as a 
function of the parameter x are given by: 

m(x) = m + x(m p h ys . — fh) (37) 

where 777 is the meson mass in the symmetry limit. 

Two different values of rTi were used, for the light mesons (the ones belonging to 
the octet) it was set to 430 MeV and for the heavy ones, 1900 MeV. 

Also the correction parameters were changed along with x, although they just 
violate SU (4) symmetry: 

7(x) = 1 + Xi^phys. - 1) (38) 

Similar functions are constructed for ^3 and ^5. 

All scalar resonances in the same multiplet have the same mass once SU{3) is 
restored while its breaking splits the masses of the different isospin multiplets. So, 
when written in the SU{3) basis, the non-diagonal elements of the matrix V (the ones 
which represent mixing between different SU (3) multiplets) are always proportional 
to 777^ — m 2 K . Figures [3] and [4] show the results for the C=0 and C=l sectors, varying 
x from to 1 in steps of 0.2. We should note that some resonances, for example, the 
k and the Dq(2400) appear as cusps for small x. This happens because thresholds 
appear during the symmetry breaking procedure. 

Table [2] displays the experimental situation of the scalar resonances, Table [3] 
shows the results in the open charm sector for the problem solved with the chiral 
Lagrangian and Table H] shows the pole positions found within the phenomenological 
model developed in this work. 

In the following the results for each sector will be discussed separately. 

5.2 C=0, S=0, 1=0 

Our model successfully generates poles which can be associated with the known light 
scalar resonances. In this sector, in the low energy region, two poles can be found in 
the T-matrix, one corresponding to the /o, but with a lower mass than one expects 
and another one for the a. It is actually possible to adjust the mass of the /o pole 
in our model by increasing the ol parameter, but two prices are paid: first the ao 
pole in the S = 0, I = 1 sector disappears for much bigger ai and also the width of 
a more massive /o decreases. 
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Figure 3: Light Scalars. The octet starts from E=851.76 MeV and the singlet from 663.13 
MeV. 
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Figure 4: Heavy Scalars. The antitriplet starts from 2327.96 MeV and the sextet from 
(2394.87-i 219.33)MeV 
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Resonance ID 


C 


s 


I 


Mass (MeV) 


T (MeV) 


fo 











980±10 


40-100 


a 











400-1200 


250-500 


a 








1 


984.7±1.2 


50-100 


K 





1 


1 

2 


841±30i^ 


618±90^ 4 


D; (2317) 


1 


1 





2317.3±0.4±0.8 


< 4.6 


D*(2A00) 


1 





1 

2 


2403±14±35 


283±24±34 










2352±50 


261±50 



Table 2: Data from [36 



Resonance ID 


C 


s 


I 


RE(v/i) (MeV) 


IM(^) (MeV) 


D: (2317) 


1 


1 





2315.41 





D*(2400) 


1 





1 

4- 


2147.65 


-107.29 


(?) 


1 





2 


No pole 




(?) 


1 


1 


1 


2427.70 


-248.40 


(?) 


1 


-1 





2410.26 


-193.80 



Table 3: Pole positions for the chiral Lagrangian 



Two more poles can be expected in this sector, one from the octet and the other 
one from the singlet, coming from the scattering of the heavy mesons. For x=0 
both poles appear, the singlet always very narrow, because its coupling to the light 
channels is very suppressed, and the octet with a much bigger width. The octet 
state disappears before x = 1.0. One should notice that the width found for this new 
heavy resonance is very small. This happens because, as mentioned (see Table [5]), 
the couplings to the light channels are very suppressed and the other possible decay 
channel is an octet formed by r/ c with a light meson which violates SU(3). Table [5] 
shows the absolute value of the residues for the resonances in this sector. 

Figure [5] shows the absolute value of the square of the transition matrix for this 
sector, as an illustration. 

5.3 C=0, S=0, 1=1 

In this sector the model successfully generates the ao resonance. Both the mass 
and the width found for it in the model agree very well with experimental values. 
Note, however that this sector was actually used to fit ol, but fitting just this one 
parameter, both the width and the mass for the ao are in good agreement with 
experiment. As mentioned if we used the pole position of the /o resonance to adjust 
the parameter at we would lose the ao pole. This relative instability of the ao 
resonance with respect to the parameters of the theory is not new, it also occurs 
when using the inverse amplitude method for unitarization and the potential of the 
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Resonance ID 


C 


s 


I 


RE(v^) (MeV) 


JM(y/s) (MeV) 


fo 











916.83 


-19.28 


a 











610.33 


-144.10 


(?) 











3718.87 


-0.11 


a 








1 


987.45 


-38.40 


K 





1 


1 

2 


840.91 


-127.48 


D: (2317) 


1 


1 





2317.32 





L>*(2400) 


1 





1 

4- 


2129.22 


-156.95 


(?) 


1 





2 


2690.50 


-439.09 


(?) 


1 


1 


1 


2700.15 


-456.88 


(?) 


1 


-1 





2686.13 


-430.87 



Table 4: Pole positions for the phenomenological model 



lowest order chiral Lagrangian where the ao appears as a cusp and not a pole. The 
pole is, however, regained when the information of the second order Lagrangian is 
used as input in the potential [23] . 

Table [6] shows its couplingqj to the different channels. 

In the heavy sector again the pole for the octet just appears for small values of x. 

5.4 C=0, S=l, 1=1 

The pole generated here should be identified with the n resonance. This resonance, 
however, is a very broad one and although there is debate on the existence of this 
resonance, recent experiments have come to support it |42|-|47], 
Again there are no heavy resonances for x=l. 

5.5 C=l. S=l, 1=0 

The D*q(2317) is reproduced in this work as a mixed bound state of \DK > and 
\D s r] >. Experimentally the observed decay channel for this resonance is D s tt which 
is not allowed in the model because it is an isospin violating process. However if 
one considers isospin violation by solving the Bethe-Salpeter equation for charge 
eigenstates instead of isospin ones and considering the real masses of the mesons, 
including the differences between different ^3 components, one gets a very narrow 
width of less then a keV for this resonance. Another possible source of contribution 
is to consider tj — tt mixing by means of which in [17] one gets a width of the order of 
a few keV. The width of this resonance is given as an upper bound of about 4 MeV 
in [36] , 

The couplings of this pole to the channels is shown in Table [8] for both models 
considered in this work. 

2 because of the identical particles, the irir channel in I = 1 just contributes to odd parity partial waves 
(indeed, the amplitude has a p-wave structure t — u). 
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Channel 


f res (GeV) 


a res (GeV) 


Heavy Singlet res (GeV) 


7T7T 


1.96 


4.23 


0.21 


KK 


3.82 


1.28 


0.03 


f]r] 


4.47 


0.47 


0.00 


DD 


0.71 


4.08 


10.41 


D S D S 


3.73 


0.49 


6.73 


VVc 


2.07 


1.04 


0.29 



Table 5: Residues for the poles in the C=0, S=0, 1=0 sector 



5.6 C=l, S=0, I=| 

Two poles are found here, one is the antitriplet companion of the -D* (2317), also 
experimentally known and to be identified as -Dq(2400). Although the antitriplet pole 
generated by the model in this sector has a width in agreement with the experimental 
value, the model fails in predicting its mass by around 150 MeV, which might not be 
too serious considering that the experimental width is around 300 MeV. 

Additionally another state is generated, belonging to the sextet. Here the two 
models differ from each other. In the chiral model this resonance has a smaller mass 
and width, but disappears as x reaches 1 because of thresholds effects, while the pole 
is predicted by our model around 2700 MeV but with a huge width that makes it 
irrelevant from the experimental point of view. 

Residues for the -Dq(2400) pole are in Table O for both models. 

5.7 C=l. S=l, 1=1 and C=l, S=-l, 1=0 

The other two states belonging to the sextet are to be found in these sectors. However 
they differ in mass and width from one model to the other. While with the chiral 
Lagrangian these poles have mass around 2400 MeV and width about 0.5 GeV, within 
our model their mass is 300 MeV larger and the huge width of the order of 1 GeV 
would make these poles irrelevant from the experimental point of view. 

5.8 Comparison With Other Works 

The light scalar resonances reproduced in this work have been thoroughly investigated 
in more sophisticated approaches and with higher orders of the chiral Lagrangian 
[22] . [23] . [24] . [28] . [29] . In our study of the hidden charm states we have now used 
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HI 2 




Figure 5: TT^ for 7r7r-channel in C=0, S=0, 1=0 sector 



coupled channels involving light and heavy pseudoscalar mesons and we find actually 
a negligible mixing of the two sectors. 

The open charm sector has been studied by Kolomeitsev [15] and Guo [T7] in 
a very similar framework but with different Lagrangians from ours; both have used 
the same Lagrangian, and very similar parameters. In [16] higher order chiral La- 
grangians are used in this sector. The Lagrangian in these works neglects exchange 
of heavy vector mesons while the present work includes it although suppressed in a 
proper way. The second term of the Lagrangian in eq. (]22h can be identified with the 
lowest order chiral Lagrangian used in [15 J and [T7] except that in the present work 
this term of the Lagrangian is a factor | smaller. Another difference between this 
present work and previous ones is the meson decay constant, /. In previous works 
it was always set to the pion decay constant, while in the present one, inspired by 
experimental measurements and lattice calculations we use a different value for the 
decay constant of the charmed mesons. 

In the S=l, 1=0 sector the results of all works coincide and the -D* (2317) is 
well reproduced. Its anti-triplet companion, the Z)q(2400), is also well reproduced in 
the S=0, 1=^ sector. However, in this sector the present work differs from previous 
ones: while within our model, the sextet state is extremely broad, in the works of 
Kolomeitsev and Guo a narrow state is predicted in this sector. The chiral Lagrangian 
we used seams to give an intermediate situation between our work and these previous 
ones, it generates for the sextet states a broad resonance although not as broad as in 
our model. The huge width of these resonances within our model is also a consequence 
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Channel 


a res (GeV) 


7T7T 





KK 


3.84 


r/Ti 


2.65 


DD 


3.64 




1.63 



Table 6: Residues for the a pole 



Channel 


k res (GeV) 


Kn 


3.81 


Kr] 


2.00 


D S D 


4.10 


K V c 


2.25 



Table 7: Residues for the k pole 



of its much bigger mass which causes a much bigger phase space for decaying into 
the open channels. 

Another novelty in the present work is the study of the hidden charm sector. Here 
we mixed light with heavy pseudoscalar pairs and concluded that there was barely 
any mixing of the heavy and light sectors. This result supports the findings for the 
light scalars, using only light pseudoscalar mesons as building blocks. On the other 
hand we find a heavy scalar with mass around 3.7 GeV corresponding mostly to a 
DD state, very narrow (see Table H]), which is surprising in view of the large phase 
space available for decay into pairs of pseudoscalars. The dynamics which prevents 
the mixing of the heavy and light sectors is responsible for this very small width. 

We should also note that with a different formalism using the Schrodinger equa- 
tion with one vector-meson exchange potential, DD states also appear for some 
choices of a cut off parameter in [48]. 
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Channel 


Chiral 
model res (GeV) 


Phenom. 
model res (GeV) 


DK 


10.21 


10.36 




6.40 


6.00 


D S T} C 


0.48 


1.52 


Table 8: Residues for the D; o (2317) pole 


Channel 


Chiral 
model res (GeV) 


Phenom. 
model res (GeV) 


Dtt 


8.91 


10.87 


Drj 


1.36 


3.77 


D S K 


5.71 


8.52 


Dr)c 


3.20 


2.14 



Table 9: Residues for the D*(2A00) pole. 



6 Summary 

We studied the dynamical generation of resonances in a unitarized coupled channel 
framework. We constructed a Lagrangian based on SU(4:) flavor symmetry and after 
decomposing the field of pseudoscalar mesons in this Lagrangian into its SU (3) com- 
ponents, we were able to identify terms mediated by exchange of heavy vector mesons 
and thus suppress these, hence breaking the SU (4) structure of the Lagrangian. The 
results were also compared with previous works based on chiral theory and heavy 
quark symmetry and with results obtained from a chiral Lagrangian considering fla- 
vor symmetry breaking effects. 

The amplitudes calculated from this Lagrangian, written in a SU(3) basis, show 
in which sectors the interaction is attractive so that it might generate resonances. 
Within the framework developed in the present work an SU (3) octet and a singlet of 
scalar mesons appear in the light sector. These resonances can be identified with the 
light scalar mesons, a, fo, ao and k, which have been thoroughly investigated before, 
but which also show up, practically undisturbed, in the enlarged basis of coupled 
channels used in the present work. 

In the heavy sector an anti-triplet is generated leading to two states which can 
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be identified with the controversial D* (2317) and Dg(2400) states, though the mass 
generated for this second one is somewhat lower than the experimental one. Thus, in 
the framework developed here, these scalar states should be interpreted as bound and 
quasi-bound states in coupled channels: The Z?* (2317) being mainly a DK bound 
state with no decay, except for a tiny one when allowing isospin violation and the 
-Dq(2400) a Dir resonance. 

Also a very broad sextet is generated in the heavy sector, but these states are 
extremely broad, making them irrelevant from the experimental point of view. One 
should note, however that these broad states contrast with states generated in pre- 
vious works [15] . [17] where narrow structures are found with the same quantum 
numbers. The Lagrangian in these previous works neglects the exchange of heavy 
vector mesons and uses a much stronger coupling. 

Also a heavy singlet appears as a pole in the T-matrix. This singlet comes from 
the attraction generated in the 3 (g> 3 sector and its structure is mainly a DD quasi- 
bound state. 
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A Amplitudes 



This appendix shows the amplitudes obtained from the Lagrangian in eq. (I29p . In 
the column of the states, the following momenta assignments should be taken into 
account: where it reads M 1 M 2 -> M 3 M 4 it means M 1 (p)M 2 (k) -» M 3 (p')M 4 (k') and 
the Mandelstam variables are defined as follows: 



a = (p + k) 2 = (p' + k'f (39) 
t= (p-p'f= (k-k'f (40) 
u= (p-k') 2 = {k-p'f (41) 

When inserting these amplitudes (or transformed to isospin or SU(3) basis) in 
the BS-equation, one should be careful to divide the amplitude by 4g each time the 
initial or the final state contains a pair of identical particles, (unitary normalization) 
in order to ensure closure of the intermediate states. The extra normalization for the 
external lines must be kept in mind but does not matter for the pole search. The 
factors 7, ip3 and ^5 are defined in eqs. d26j) - (J2HJ) . 



A.l C=2. S=2 



States 


Amplitude 


D+D+ -> D + D + 

s s s s 


-^pi-M^s ~t-u) + 2m 2 D + 2m 2 K - 2m 2 n ) 


A.2 C=2. S=l 


States 


Amplitude 




-Wp{-^(s -u)-{s-t) + 2m 2 D + m 2 K - ml) 


A.3 C=2. S=0 




States 


Amplitude 




D + D » D + D u 


-u)-(s-t)+ 2m 2 D ) 



A.4 C=l. S=2 



States 


Amplitude 




- «) - l{s -t)+m 2 D + 2m 2 K - m 2 ) 
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A.5 C=l, S=l 



States 


Amplitude 


XV iy " 11 1 J 

_». if°.D+ 


epwl' "J t ^ u>) -r iii D -f J 
i,*- (n/(f — 7/^ -t- f s — 7/^ -t- m -1- m?. 

6v 4 /2 ( (« «) 7(s t)+m^ + m^) 
_ 6^^ - t) - (3 + 7 )(s - u) - rr& - 3m|- + 2m;) 




-6p(7(* - U ) + ( S - U ) + m D + 

6v 4 /2 (( s ^) + 7(s t) m 2 D m\) 


-^r,D+ 




r)D+ - t/Z>+ 


-97r(7(-a + 2t - n) + 2m^ + - 4m£) 


r/ c L>+ -> 


-T872-(4t(-s + 24 - it) + llm£ + 3m 2 K - 7m'l) 
-^kp^~ s + 2t ~ u )~ 2m D + O 
~^kp^~ s + 2t ~ u )~ 2m l + ml) 

-^kw^- 2 <~ s + 2t ~ u ) + 2m D ~ ml) 


A.6 C=l, S=0 


States 


Amplitude 


-► 7r _ .D + 

- K-Dt 


-ji]?z{l(-s + 2t-u) + 2m 2 D + 2ml) 

-i^(7(-s + 2t-u) + 2m 2 D + 2 ml) 
6v4f2 u ) + s u + m D + m\) 


TT-D+ -► 7T-D+ 


-6F(7(* - u) + s - u + m 2 D + ml) 
-^kp^~ s + 2t ~ u ) + 2m D + 2 ml) 
-■&fi{l{t - u) + s-u + m? D + m\) 




-36^(7(-s + 2t-u) + 2m 2 D + 2m 2 n ) 
~^kp^ s -t) + (3 + 1 )(s-u)-m 2 D - 3m 2 K + 2m;) 


X _ £>+ K~D+ 


~6f(7(* - + s - n + m|, + 2m? K - m 2 ) 


Vc D u -► t? cj D u 

— > 7T°Z) 
-► 7T-D+ 
-► f?L>° 
- ^~^ + 


- T gj 7 (47(-s + 2t - u) + llm^ - 4m^) 
" W^(" s + 2t ~ u )~ 2m D + 

" W (27( " S + 2t " n ) " 2m ?) + m ') 
-i^(27(-s + 2t - «) - 2m 2 D + ml) 

-67sp( 2 ^- s + 2t ~ u )~ 2m l + ml) 
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A.7 C=l. S=-l 



States 


Amplitude 


K-D+ 

— > K°D° 


- u) - 7 (s - t) + rr& + m^) 







A.8 C=0, S=l 



States 


Amplitudes 


D+D" D+D- 


-6^(* ~ n + ^(s - u) + 2m|, + m|- - m£) 


-► 7T-K+ 


- i^r(-s + 24 - u + 2m^ + 2m£) 
2v 4 /2 ( * + ^) 
"TiTsT^ 5 " 2 * + u ) + 2 < " 2m ') 




-gj?(s + t - 2u + + m 2 n ) 
-^kp(-^ s ~ 2t + u )- 2m K + 2m*) 


r/K u -► ^ 


- ijjr(-3(a - 2t + «) + 6m^ - 2m£) 


rfcif u rfcif u 


1 2 


D+D" - vr u ^ u 


6v 4 /2 (* u 7(3 *) ™D m l) 

-6^(-(t -u) + 7(5 - t) + mfj + m^) 
_ 6^« 3 + 7)(« - *) - 7(a - u) - m£ - 3m| + 2m£) 


D+D" -► ?7 CJ FC U 


~6tM 2 ^ 2s -t-u)-2ml + mi) 




6v 4 /2 ( mi) 
6v | /2 (2< + mj ) 
6 ^ f2 ( 2m| + m2) 
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A.9 C=0, S=0 



States 


Amplitude 


D+D- - DfD- 
-► 


-37^3(5 + * - 2u) + 2m7j + 2m^ - 2m£) 
-■Az{ip5(t -u) + s-u + 2m 2 D + m\ - m\) 
-^z(4>5(t -u) + s - u + 2m? D + m 2 K - ml) 


i j ij — > ij i j 


L filial a _1_ + Oni\ _1_ 9tti^ ^ 

— 1^3 (s -\- 1 — zu ) -\- zm D ) 
--^{^{t -u)+s-u + 2m 2 D ) 




-^Ms + t-2u)+2mi ) 


K + K~ -► K+if- 

— > 7T~^7T 

-ttV 
— > 7r°r7 
— ► 7717 


-jp{s + t-2u + 2m 2 K ) 
--± 7 {s + t-2u + 2m 2 K ) 

1 ( e 1 + On J- J- m^A 
— qJ^V* • ^ — ^ W ' K ' ''''nJ 

-j^p(2s -t-u + 2m\ + 2m 2 ) 
-JtTspW 28 -t-u)-2m 2 K + 2ml) 
-j^(3(2s -t-u)+ 6m 2 K - 2ml) 


K U K U — > K U K U 

— > 7T + 7T~ 

-ttV 

n 

-> 7T°v7 
— ► 7777 


~3p\ s + t — 2u + 2m K ) 
-6p( s ~2t + u + m\ + m 2 ) 
-j5ji{2s -t-u + 2m\ + 2m 2 ) 
"i^(-3(2s - t - u) + 2rr& - 2ml) 
-1^(3(25 - t - u) + 6m 2 K - 2ml) 


7T + 7T~ — > 7T + 7T~ 

— > 7T°77 
— ► 7777 


-^p(s + t-2u + 2mi) 
-jp(2s-t-u + ml) 


7r u vr u 7T U 7r U 
— > 7T°?7 

— ► 7777 


—prrii 
1 2 


7T U 77 — > 7T U T7 
-> 7777 




7777 — > 7777 


-^(16m 2 -7m 2 n ) 


v7 c 7T° -> ?7 c 7r u 


1 2 






-^K°K° 
— > 7r + 7r _ 
-ttV 
— > 7r°r7 
-► 7777 


- 6p(* - « + 7(s - «) + m^, + 2m^ - rri*) 
- u + 7(s - ti) + mjj, + 2m^ - m^) 

-gi7(7(2s -t-u) + 2m 2 D + 6m|- - Am 2 ) 
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States 



Amplitude 



D+D~ 



K + K- 

-> TT + TT" 

^7r vr° 

— > 7T ?7 
-»• 7/7/ 



" 67*(-(* - U ) + 7(S - *) + ™D + m \) 

-Gj?(t -u + j{s-u)+m 2 D + m 2 ) 
-^(7(25 - i - u) + 2m 2 D + 2 "4) 
_^_(_ 7 ( 2s _ t - u) - 2m% ~ 1ml) 

"36^(7(25 



t-u) + 2m 2 D + 2m l) 



jz(-(t -u)+ 7(s -t) + m'l + m\) 



K + K- 

-> 7T + 7T~ 

-►A 

-> 7/7/ 



6/ 



6/ 



i 7 ( s - *) + m 2j + m 2) 



^(7(25 -t-u) + 2m 2 D + 2 "4) 



_^_( 7 ( 2s _ t - u) + 2m 2 D + 2<) 
-3^(7(2 S -t-7i) + 2m| ) + 2m^) 



^7 



77 C 7T" 

-»■ »7c»7 



-> 7/ c ?7 



■ _^_(_2 7 (2 S -t-7i) + 27721, 



7jj + 2m 2 D — m. 
u) — 2m D + m\ 



D + D~ 



-»■ ?7c»7 



ev^ 2 



(-2 7 (2s -t 
(27(2s - i - 



2ml) + ml) 
- 2ml + ""4) 



18v^2/- 



(2 7 (2s - t 

:(2 7 (2 S -/ 



1 ( 2ttj 2 - — 2 



1 2 

' 3^/2/2 m n 



7/ c 7T^ 



K+K- 

-> 7T + 7T~ 

— > 7r 7/ 
— > 7/77 



m 2 ) 



■ ^(-2m^ + m 2 ) 



7/ c 7/ — > K + K 
— ► 7T + 7T~ 

-ttV 

— ► 7T ?7 
— ► 7/7/ 



1 2 

' 372/2 ^ 
1 2 

" 372/2 ^ 



~97f?2- ( ~ 8m ^ + 5m ^ 



B Isospin and SU(3) basis 



The following phases are taken for the meson assignments of the 15-plet: 

\D + > \ _ ( \K+> \ 

^ -\d°> y |K > \- [\k°> y 



\D S > =\D+>, \D> V 



7T >l = 



|7T U > 



|r? > = h >, \Vc >o= |% >, l-K" >l 



-\K~ > 
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\D° > \ 

\D>i=[ ' _ ^ and \D S > = \D S > 



2 



jzr > 7 

In the following we will list for the sectors where a SU(3) decomposition is not 
trivial, the isospin and SU(3) states used to transform the amplitudes from a charge 
basis to isospin and then from isospin into a SU(3) basis. For reviews on phase 
conventions and isoscalar factors of the SU{3) Clebsch-Gordan coefficients one can 
refer to OS], EDI. 



B.l 3® 3 (C=2) 

\d s d s 




\D+D° > 
\D°D+ > 




Figure 6: 3 <g> 3 = 3 © 6. 



|6,2,0 >H= \D S D S > 
\6,1,±>=±(\DD S >+\D S D>$ 

|6, 0, 1 >= \DD > 

\3,1,0 >=^(\DD S >-\D S D>) 

|3,0, \ >= \DD > 




C=l 



\K°D+ > 



KD > \ = j_ ( -1 -1 W \K+D° > 
KD> 1 J V2 \ -1 1 J y \K°D+ > 

\ V D S > = \ V D+ > 
\ttD s >!= \ir~D+ > 



SU(3) states are represented as \Irrep, S,I > 

from now on the label for the isospin of the states will be omitted for the SU (3) states. 
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[\nD> 3 J ^ ){\«-D + >) 
\r)D >i= -\rjD° > 
\KD S >i= -\K~D+ > 

\KD > \ = j_ ( 1 -lW \K-D+ > \ 

\KD >i j V2 y -l -l J y \K°D° > J 




Figure 7: 8 <g> 3 = 6 © 15 © 3. 

|f5,2,± >= \KD S > 

|15,1,1 >= -fe{\KD > -\ttD s >) 

115, 1,0 >= -%VD S > +\\KD > 
|f5,0, § >= \itD > 

1 15,0, \ >= \\itD > + \\qD > -yJ\\KD s > 
1 15, —1, 1 >= \KD > 
\6,1,1>=±(\KD>+\ttD s >) 

|6,0, \ >= j\\*D > -y/%\r,D > -\\KD S > 

|6, —1,0 >= \KD > 

|3, 1,0 >= \\r)D s > +%KD > 

|3,0, i >= -\\-kD > > -J\\KD S > 

B.3 3® 3 

\D S D >i= |D+£>" > 
|D a J fl > = \DfD- > 
( \DD> \ = ^( 1 
^ \DD >i y v/2 ^ l 

|8,l,i >= |L> S D> 
|8,0,1 >= |L>L> > 

|8,0,0 >= J\\D S D S > 

|8, -1, | >= \D S D> 



c=o 



1 \ / \D+D- > 
-1 J I |L>°L> > 



-73^ > 
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Figure 8: 3 © 3 = 8 © 1. 



|1,0,0>= ^\D S D S >+Jl\DD> 

B.4 8® 8 C=0 

(\« K >\\ = ( 7a ->/f V l*°*°> 
|r?K >i= |r/K° > 



I™ >o ^ _ JL ( - 1 \ K+K ~ > \ 

\KK >! J V2\-l 1 J \ \K°K° > J 




Figure 9: 8 <g> 8 = 1 © 8 © 8 © 10 © 10 © 27 
|27,2,1 >= |KK > 
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27,1,1 
27,1 
27,0 
27,0 

27,0 

10,1 
10,0 
10,1 
10,0 

8 5 ,1 
8 5 ,0 
85,0 

8 A , 1 
8 A , 
8a, 



>= 

h>= 

2 >= 
1 >= 

>= 

h>- 

1 >= 

1 > = 
1 > = 

>= 

1 >= 

>= 



^(|Kvr > +\ttK >) 

^(I^tt > +\ttK > -3\Krj > -3\r]K >) 

1 7T7T > 

> +\KK >) + -^(Ivrr? > +\ W >) 

> ~\ RK » + ^k> "iTsolw > 
i(|K7r > — Itt^T > +|Kr? > -|??^ >) 

> —\KK > -\irir >) + \{\irr} > -\r]ir >) 

^(-lETvr > +\ttK >) 

4=(-|i<r^ > +\KK > +\tttt >) + \(\ttt] > -|r/7r >) 
: ^{-3\Kn > -3\itK > -\Kr] > -\rjK >) 
^{-\KK > —\KK >) + ^(Ivrr? > +\rpr >) 

^o(-|™ > "I™ » + ^ > +7SIW > 
: 1(-\K-k > +\itK > +\Kr] > -\r]K >) 

■jg(\KK > —\KK > +2|vr7T >) 
-^(-\KK > +\KK >) 



1,0,0 >= k{-\KK > —\KK >) 



2 75l 7r7r>+ 272lw > 



C Isolating The Contribution from £3 

The £3 Lagrangian in eq. (|20p has two terms, the first one contains just hadronic 
currents where the initial and final state have the same electric charge and, therefore, 
exchange neutral vector mesons only. The other term has contributions from both 
charged and neutral vector mesons; from this second term first one should isolate the 
contribution from neutral vector mesons. Let us add and subtract the appropriate 
term, 

Tr ( J 33/^33) Tr ( y -h3vJ^-{ y J D s D a v J D s D s +J D+D-v J D+D- +J D°D°^ 

+ ( J D s D s n J D s D s + J D+D-v J D+D- + J D°D °ll J D°D<>) ( 42 ) 

such that now the sum of the first two terms in eq. (|42p has no contribution from 
heavy vector meson which is now in the third term alone. 

The second term of Lagrangian £3 in eq. (|20|) will then be expanded in order to 
identify terms where equal hadronic currents are connected and terms where different 
ones are connected: 

J 33^3 = 2 {Jd s D s li{Jd+D- + ^D°D°) + J D+ D~ ^J^o fro) + 
J D a D af j, J D s D a + Jd+D-^Jd+D- + ^£> -D°M^D°S0 

Now terms with the product of equal neutral hadronic currents are to be multi- 
plied by the correction ^3 and terms connecting different ones by ^5, given in eqs. 
(gZD and pgj. As a result: 
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Jd+d-h Jd+d- + Jd°d°h Jd°d°) 
One can work it out: 

(2^3 - 1) ( J DsDs^ J D s p s + J D+D-n J D+D- + ^CSV J D°5°) = 

~>p5 

Tr (^33^^35) + ^[JdsDs^DsDs + J D+D-v J D+D- + JdODO^DODO + 
2Jd s D s h(Jd+D- + ^D°S°) ~*~ 2 Jd+ D-^i J popo) = Tr (^33^^33) + ^5^33/2^3 

And this is the simple form we write down in eq. (|29l) . Yet, in the amplitudes we 
use the ^3 and V>5 factors. 
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